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Abstract: We calculate the heavy quarkonium energy levels and decay widths in a quark- 
gluon plasma, whose temperature T and screening mass m£, satisfy the hierarchy mas ^ 
T 3> mOg >> ruD {rn being the heavy-quark mass), at order ma^. We first sequentially 
integrate out the scales m, mag and T, and, next, we carry out the calculations in the 
resulting effective theory using techniques of integration by regions. A collinear region is 
identified, which contributes at this order. We also discuss the implications of our results 
concerning heavy quarkonium suppression in heavy ion collisions. 
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1. Introduction 

Heavy quarkonium has been suggested since long time as a thermometer for the medium 
that forms at the core of heavy-ion colhsion experiments ^. The early arguments were 
based on the naive expectation that above the deconfinement temperature the confining 
part of the quark-antiquark potential vanishes and the Coulomb part turns into a Yukawa 
potential due to screening. Since the Yukawa potential supports a finite number of bound 
states depending on the screening (Debye) mass, and the latter is linear in the temperature, 
it is then clear that the relative fraction of the different heavy quarkonium states observed 
will depend on the temperature of the medium. In addition, the electromagnetic decays 
of these states provide a clean experimental signature. The gross picture above appears to 
be supported by experiments [^. 

In the last few years, significant progress has been made in deriving the quark-antiquark 
potential on a rigorous basis. A model independent study of the real-time static potential 
was initiated for large temperatures (T S> l/r > mu) in [y, Q, |5|, g] and its implications 
for a QED and QCD plasma studied. For a wider range of temperatures, including lower 
temperatures, an effective field theory (EFT) study of non-relativistic bound states in a 
plasma was initiated in for QED and in ||8[ for QCD in the static limit. The potential 
obtained in this way differs in many respects from the most commonly used phenomeno- 
logical potentials (for some reviews see |^, |l^ ) . Most remarkably, it develops an imaginary 
part. At least two mechanisms have been identified that are responsible for the appearance 
of a thermal width: the Landau-damping phenomenon Q and the quark-antiquark colour- 
singlet break up Q. In particular, it has been pointed out that quarkonium dissociation 
due to the former rather than screening may be the dominant mechanism at the origin 
of heavy quarkonium dissociation in a medium j^, [ill. These developments motivate us 
to revisit the physics of heavy quarkonium states in a thermal bath in a more systematic 
way. We shall focus here on temperatures for which ttT is smaller than the typical momen- 
tum transfer in the bound states: such temperatures are those reachable at present days 
cohiders i|lO|. 

Heavy quarkonium in a medium is characterized both by the scales typical of a non- 
relativistic bound state and by the thermal scales. The non-relativistic scales are the 
inverse of the typical radius of the system 1/ao and its typical binding energy E. The 
thermal scales are the temperature (or multiple of vrT) and the electric screening mass 
rriD, among other lower energy scales, which are not relevant to our discussion. In the 
weak-coupling regime, which we will assume throughout this work, these scales may be 
expressed in terms of the strong coupling constant g <^ 1, the heavy quark mass m, and 
the temperature T: m^ ~ gT, 1/ao ~ mus and E ~ mOg, where Og = 5^/(47r). Non- 
relativistic scales and thermal scales are hierarchically ordered. This allows to investigate 
the quarkonium properties in a medium using the same systematic framework provided by 
non-relativistic effective field theories at zero temperature |12]. 



In this work, we aim at studying heavy quarkonium at finite temperature including the 
contribution induced by a large but finite quark mass, in this way merging and completing 
the findings of Refs. ||^, ^]. We will adopt the same real-time EFT framework of 0, |8| and 



assume for definitiveness the following hierarchy between the thermodynamical and the 
non-relativistic scales: 

m » mus ^ T » ma^ S> mn. (1-1) 

With this choice, the thermal bath affects the Coulombic bound state as a small perturba- 
tion, yet modifying the Coulomb potential. We remark that this temperature is below the 
melting temperature, which is of order maj R. Moreover, this may indeed correspond 
to the situation of interest in present day colliders. As a consequence of ( |1.1| ), in the weak- 
coupling regime, we have that mg^ S> T ^ fng'^, corresponding to mg^ ^ m/) ^ mg^ . We 
furthermore assume that Aqcd, the QCD scale, is smaller than mo (although results that 
do not involve a weak-coupling expansion at the scale mo , which are all the results of the 
paper before Sec. |6|, are valid also for mo ~ Aqcd)- A number of different inequalities has 
been addressed in the Abelian case in [13|. 

We will concentrate on the energy levels and decay widths and we will determine 
how they get modified in a thermal bath whose temperature is such that it satisfies the 
conditions ( |1.1| ). In order to be definite, we will further assume {m£)/E)^ <^ g, in this way 
keeping small the number of required corrections suppressed by powers of mo/ E, and we 
will evaluate the spectrum with an accuracy of order ma^ . 

The strongest limitation for the practical application of our final results to actual 
bottomonium and charmonium systems comes from the fact that we use perturbation 
theory at the ultrasoft scale mag . Still, we expect them to be relevant for the ground states 
of bottomonium and, to a lesser extent, charmonium. Some intermediate expressions, for 
which perturbation theory is only used at the scale T ^ mal may have a wider range 
of applicability. We also assume a vanishing charm quark mass in the bottomonium case 



(effects of a non- vanishing mass are discussed in [13|). 

The paper is organized in the following way. In Sec. |2|, we briefiy review the Feynman 
rules of QCD at finite temperature in the real-time formalism. In Sec. y, we set up the 
effective field theory that follows from QCD by integrating out the scales m and mag in the 
heavy quark-antiquark sector. In Sec. ^, we calculate the contributions to the spectrum 
coming from the scale T, in Sec. ^, those coming from the scale E and, finally, in Sec. |6|, 
those coming from the scale tud- In Sec. ^, we summarize our results giving the thermal 
energy shifts and widths up to order mOg. 

2. QCD at finite temperature in the real-time formalism 

In this section, we review the Feynman rules of QCD with static quarks at finite tempera- 
ture in the real-time formalism. 

Real-time expectation values depend on how the contour of the time integration in the 
partition function is deformed to include real times. In the paper, we adopt a contour that 
goes from an initial time tj to a real final time if, from if to tf — iO+ , from tf — zO^ to t\ — iO^ 
and from t; — iO'^ to t\ — i/T. The propagators will be given with this conventional choice of 
contour. Since the contour has two lines moving along the real time axis, degrees of freedom 
double in real time and propagators are represented by 2 x 2 matrices. Furthermore we 



define 



nB{ko) 



1 



(2.1) 



gfco/T _ 1 ' 

which is the Bose-Einstein distribution. 

The non-relativistic propagator of an unthermahsed quark-antiquark pair interacting 
through a potential V{r) reads (see S, we have added here the kinetic energy) 



S{ko,k) 



ko — k'^/m — V{r) + irj 
2TT6{ko - k'^/m - V) 







(2.2) 



ko — k'^/m — V{r) 



trjj 



The expression for V{r) depends on whether the quark-antiquark pair is in a color singlet 
or in a color octet configuration. In the last case, an identity matrix in color space must 
be understood, and, in either case, an identity matrix in spin space is implicit. Since 
the \S{kQ,k)\i2 component vanishes, the quark-antiquark fields of type "2" never enter in 
any amplitude of physical fields, i.e. fields of type "1". As a consequence, the fields "2" 
decouple and may be ignored when considering physical amplitudes. 

Throughout this paper we adopt the Coulomb gauge for our calculations. The free 
gluon propagator reads 114 



D^^ik) 
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+27T5ikl-k^)nB{\ko\) 




(2.4) 



where k stands for the modulus of the three momentum fc* . 

Due to the decoupling of the heavy quarks of type "2" , in the following we will need 
only the "11" component of the heavy quark-antiquark propagators. At the order we are 
calculating, this is also the case for the gluon propagators, which will be needed up to one 
loop. All our equations will refer to this component, unless explicitly stated otherwise. In 
particular, we recall that at equilibrium the "11" component of the gluon propagator can 
be written in terms of the retarded (R) and advanced (A) propagators. 



(2.5) 
(2.6) 

(2.7) 



as 



D%{k) = - J d'^xe^^'^''^^^-'--^ e{-xo){[A^ix),A,iO)]), 
D^,Ako,k) + Dt,(ko,k) (\ \ , t> a 



which holds for the tree level propagator as well as for the full one. The second term 
on the right-hand side, proportional to the difference between the retarded and advanced 
propagators, is often termed the symmetric propagator. 

3. Integrating out the scales m and mas 

Out aim is to calculate the quarkonium spectrum in a thermal bath of temperature T. 



We take advantage of the hierarchy of scales (O) by constructing a hierarchy of effective 
field theories that follow from QCD by systematically integrating out the largest scales. 
The EFTs are constructed as series of operators whose matrix elements scale like the lower 
scales and that are suppressed by powers of the large scales, which have been integrated 
out. 

The first scale to be integrated out from QCD is the heavy quark mass m. In the 
matching procedure, smaller scales are expanded. Thus, at this stage, the presence of the 
thermal scales does not affect the matching of the Lagrangian, which is the Lagrangian of 
non-relativistic QCD (NRQCD) p 



The next scale to be integrated out is the scale of the inverse of the typical distance 



of the heavy quark and antiquark, which is of order mas- According to (IT), it is larger 
than the temperature. We are thus allowed to integrate out mas from NRQCD setting to 
zero all thermodynamical scales. Furthermore, under the assumption that mas ^ Aqcdj 
this integration can be carried out in perturbation theory order by order in as- The EFT 
we obtain is potential non-relativistic QCD (pNRQCD) [l^, 17|. Its Lagrangian reads 



-CpNRQCD = -^^i^"''" + X; g^ i]pq^ + J d\ Tr J S^ [iO^ - /i.] S + O^ [iD^ - /i^] O I 
+Va Tr jo^r • c/E S + S^r • 5E o| + ^Tr |o^r • 5E O + Ql'Or • g^X + ... . (3.1) 

The fields S = S" lc/\/iVc and O = 0° r"/\/Tp, are the quark-antiquark singlet and octet 
fields respectively, nj is the number of light quarks qi, Nc = 2> \s the number of colours, 
Tf = 1/2, E is the chromoelectric field {E^ = F*°) and iDqO = idoO - gAoO + OgAo. The 
trace is intended over colour and spin indices. Gluon fields depend only on the center-of- 
mass coordinate R and on time; this is achieved by a multipole expansion in the relative 



distance r. The dots in the last line of Eq. (3T) stand for higher orders in this expansion. 



The dependence on the hard and soft scales m and 1/r is encoded in the Wilson 
coefficients; Va and Vb are at leading order V^ = Vg = 1, whereas the singlet and octet 
Hamiltonians have the form 

/i.,o = ^ + yi? + ^ + ^ + ..., (3.2) 

m ' m ?7z^ 

where ttt. is the heavy quark mass, p = — iVr. The dots stand for higher-order terms 
in the expansion in 1/m, both for the kinetic terms (relativistic corrections) and for the 
potentials, as well as for terms that depend on the center of mass three momentum. 



The static potentials read 



Vi°' = -C.^, V.'»' = 4^. (3.3) 

where Cf = {N^ — l)/{2Nc) and ay^ and ay^ are series in Og and at the leading order 
ay^ = a\4 = Oig. ay^ is known up to three loops |18, ^, whereas ay^ to two loops pOJ . 



Starting from order a^, ay^ is infrared divergent. This divergence was first identified in |21| 

and analyzed in the framework of pNRQCD in |22|, where it was shown to cancel against 

an ultraviolet (UV) divergence coming from the ultrasoft degrees of freedom (the scale E') . 

The non-static potentials Vs and Vg can be read from |23, 12 1. Vg consists of a sum 



of many terms, such as a p-dependent term, terms depending on the angular momentum, on 
the heavy quark-antiquark spins and a spin-orbit term. Some of these terms, as well as V^ , 
have an infrared divergence. The leading logarithmic dependence on In(^r) accompanying 
these divergences can be read from p^, 24]. In the octet sector, the non-static potentials 



are not known beyond tree level. Fortunately, for the present analysis only the leading 
order expression in Og for the static octet potential will be needed. 

The power-counting of the pNRQCD Lagrangian ( |3.1D goes as following: the relative 
momentum p and the inverse distance 1/r have a size of O (mas), whereas the time deriva- 
tive, the gluon fields and the center-of-mass momentum P scale like the lower energy scales. 
Therefore, the largest term in the singlet potential expanded in a^ and 1/m is the Coulomb 
potential —Cpa^jr. The spectrum of the corresponding singlet Hamiltonian is given by 
the Coulomb levels 

mCpO^ 1 _ 2 
^n = 7-5 — = 2~T ' '^0 = — 7^ • i^-^) 

Subleading terms in the expansions in Og and 1/m are treated in quantum-mechanical 
perturbation theory. The corresponding shifts of the Coulomb levels have been computed 



in |23, 24, 53]. The infrared divergences mentioned above affect the spectrum at order 
mag. 

For what concerns the propagators of the singlet and octet fields in the real-time 
formalism, we have shown in Eq. ( p.2| ) that the "12" component of a quark-antiquark 
propagator in a potential V{r) vanishes and the unphysical "2" component decouples. We 
are thus allowed to drop also here the real-time formalism indices and write only the "11" 
component of the propagator. For the rest of the paper, all amplitudes will be intended 
as the "11" components of the real-time matrices unless otherwise specified. In particular, 
for what concerns the singlet propagator, we thus have 

^singlet^^) = ^ , (3.5) 

where E is the singlet energy. From Eq. (|3.2|), the singlet Hamiltonian hg reads 

hs = ^- Cf— + ■■■ , (3.6) 

m r 



Figure 1: The leading heavy quarkonium self-energy diagram. The single line is a singlet propaga- 
tor, the double line an octet propagator, the curly line a gluon and the vertices are chromoelectric 
dipoles vertices. 

where the dots stand for higher-order terms. In order to have a homogeneous power count- 
ing in the propagator, it is convenient to expand it around the leading-order Hamiltonian 
( |3.6| ) , which is of order mag . Similarly the octet propagator is 

S''^^^\E)ah= / (ite*(^-'^°)Me~^^io^o) , (3.7) 

where the octet Hamiltonian ho reads 

p^ 1 Os , ■. 

m INc r 



and the dots stand for terms smaller than mUg. The Wilson line in (^J) ^^^ ^^ expanded 
in powers of g. We will only need the leading order in such expansion, 

4. Integrating out the temperature 

In this section, we proceed to integrate out modes of energy and momentum of the order 
of the temperature T. This amounts to modifying pNRQCD into a new EFT where only 
modes with energies and momenta lower than T are dynamical. We may denote the new 
EFT with pNRQCDhtl 11 • The EFT can be used for mas > T > £;,mD no matter 
what the relation between E and ttid is. Its Lagrangian will get additional contributions 
with respect to pNRQCD. For our purposes, we are interested in the modifications to the 
singlet sector, corresponding to a thermal correction 6Vs to the singlet potential, and to the 
Yang-Mills sector, amounting to the Hard Thermal Loop (HTL) Lagrangian >Chtl |27]: 

^pNRQCDhtl = ^HTL + / d^r Tr J S^ [iOo - h, - 6V,] S + O^ [iDo -ho- 5Vo] O i 

+Tr jo^'r • ^E S + S^r • gBo\ + ^Tr lo^r • 5E O + O^Or • g^l + ... , (4.1) 

where we have set to one the matching coefficients of the dipole terms, whose quantum 
corrections are beyond the accuracy of the present paper. 

We calculate the correction SVg to the singlet potential. As in fs], 0, the leading thermal 
correction is due to the dipole vertices O^r • gE S -|- S^r • (7E O in the pNRQCD Lagrangian 
( |3.1| ). These terms induce the diagram depicted in Fig. |l|, where a colour-singlet state 



emits and reabsorbs a chromoelectric gluon through the dipole vertex and an intermediate 
colour-octet state. The amplitude reads (see |23, ^, 24] for the T = case) 



-ig^ Cf 



A~D 



D-1 



/^ 



d^k 



(27r)^ E-ho-ko + irj 



l(0) 



(0), 



k'oDi:;'iko,k) + k'D'^,'{ko,k) r\ (4.2 



where E is the energy of the singlet; we recall that this expression corresponds to the 
"11" component in the real-time formalism. Integrals over momenta are regularized in 
dimensional regularization, with D = 4: + e and // being the subtraction point. In Coulomb 
gauge, with the free propagators given in Eqs. ( |2.3| ) and ( |2.4p , the contribution of the 
longitudinal gluon vanishes in dimensional regularization, whereas that of the transverse 
gluon can be divided into a vacuum and a thermal part: 



-ig^Cp 



D 



D-1 



rV^-^ 



d^k 



kn 



{27t)^ E - ho - ko + ir) ° 



A-2 



k'^ + ir] 
2TT6{kl-k^)nBi\ko\) 



(4.3) 



the first term in the square brackets is the vacuum part and the second term is the thermal 
part. The expression depends on the scales T and E. In order to single out the contribution 
from the scale T, which comes from the momentum regions ko ^ T and /c ~ T, we recall 
that T ^ {E — ho) and expand the octet propagator as 



E-h, 



+ i 



{E - hof 



(E-ho 



E - ho - ko + iij -ko + irj {-ko + irjY (-/cq + ir}Y (~^0 + i-f]) 



+ .... (4.4) 



The contribution of the vacuum part of the propagator is scaleless for all the terms of the 
expansion and thus it vanishes. Conversely, in the thermal part, we have the Bose~Einstein 
distribution giving a scale to the integration. 

The zeroth-order term in the expansion ( |4.4| ) gives a vanishing integral |^], whereas 
the following terms contribute to the potential. The linear and the cubic terms in E — ho, 
i.e. 



-9^Cf 



D 



D-1 



r\E-hoy n 



A-D 



and 



-g^Cp 



D 



D-1 



r\E - ho 



iV/x^-^ 



d^-^k UBJk) 
(27r)^-i k ' 

d^-^k n^jk) 
(27r)^-i P ■ 



(4.5) 



(4.6) 



can be shown to contribute to the real part of the potential. Since in our counting 



behaves as mg^ ^ a^T'^Er'^ ^> mg^^ and (4.6) as a^E^r"^ ~ mg^^, further terms in the 



E/T expansion are not needed. Finally, the square term in the expansion, which would 
give an imaginary contribution to the potential, vanishes in dimensional regularization: 



g'^Cp D-2 
2 D-1 



r\E 



/io)V>^-^ 



d^-^k 
-^-yj^^kn^ik) 



1 



+ 



1 



{—k + ir])^ (k + irj)^ 



0. (4.7) 



We now evaluate the linear term defined in Eq. (^^). The integration yields 

- CpasT^r^E-Ky. (4.8) 

y 

Matching the singlet propagator in pNRQCD with the singlet propagator in pNRQCDjjrpj^ 
we obtain 

wh-. - whr. (-1 ^-.^--(^ - '•«)^') ^ - '^i^W^l^/'i'" ■ 

(4.9) 
where the left-hand part of the equality corresponds to the pNRQCD part of the matching 
and the right-hand side to the pNRQCDjjrp^^ part of the matching: 6Vs is the thermal 
correction to the singlet potential in pNRQCDjjrpj^ and SZg the thermal correction to the 
singlet normalization in pNRQCDfjrpL- Our purpose is solely the evaluation of 6Vs, which 
is necessary for the spectrum. So we rewrite E — ho as E — hg — {ho — hg), where {ho — hg) 
is given by the difference between the octet and singlet potentials: 

T/(") _ T/(") 

ho-hs = y^- -^^AV. (4.10) 

n 

AV is organized as an expansion in Og and 1/m. At the leading order, it is the difference 
between the tree-level static potentials: 

AF = ^^. (4.11) 

Higher-order terms are easily shown to contribute to the spectrum beyond our accuracy. 
Similarly, for what concerns the singlet Hamiltonian, only the leading terms displayed in 
Eq. ( |3.6D are necessary. Hence r^{E — ho)r^ simplifies into r^{E — hs)r^ — Nca^r /2; the 
second term is easily identified as contributing to 6Vs , whereas plugging the first term back 
into Eq. (|^) yields 

i rUE-hsY ^- = i ( \\r\ E - /i,], r^l + {r^, {E - /i,)}) ^ . (4.12) 

E-hs ^ ' E-hs 2E-hs^^^ si, i^x ,y "^) E-hs 

The term {r^, {E — /i^)} contributes to the normalization of the wave function. We are 
thus left with 

^ ^[[hsy]y]^r^- (4.13) 



E-hs2'' ^ J' 'E-hs 



The commutator can be easily computed from the Hamiltonian displayed in Eq. (^.6|). 
Plugging the result back into Eq. (O), we obtain from the matching condition 

^y(linear) ^ iNcCFa^T^ + —Cpa.T^ . (4.14) 

9 3m 

The first term is the contribution of AF and was first obtained in ^j. The second term is 
the contribution of the kinetic term; a similar term appears in the Abelian case of Ref . Q • 
Using first-order quantum-mechanical perturbation theory and the expectation values (r)„ ^ 
on the eigenstates of the Coulomb potential (n and I stand for the principal and angular 



momentum quantum numbers respectively, see, for instance, [p9| ]) we obtain the following 
correction to the Coulomb energy levels 



^KT = ^N,CpatT'^[3n' - l{l + 1)] + ^Cpa^T^ (4.15) 

We now move to the cubic term, as defined in Eq. ( [4.6|) . We have 

-,^C,B_lr'iE - Kfr','-- I ^'^ = ^r'(E - Kfr' . (4.16) 

where It comes from the evaluation of the integral. It reads [^ 

2 T^ 5 

/;r = - + ln^-7s + ln(4^)--, (4.17) 

where je is the Euler's gamma. The divergence of this expression is of infrared (IR) origin: 
it arises when integrating over the Bose-Einstein distribution at momenta much smaller 
than the temperature. Since we are integrating out the temperature, i.e. getting the 
contribution for /c ~ T, this divergence is an artifact of our scale separation. We identify 
two possible schemes in which the cancellation of this divergence may be interpreted. 

1. In the first scheme, the divergence is cancelled by an opposite ultraviolet divergence 
from a lower scale, in our case the binding energy. In the next section, we will 
indeed show that the thermal part of this very same diagram, when evaluated for 
loop momenta of the order of the binding energy, yields an ultraviolet divergence 
that exactly cancels the one here, whereas the vacuum part of that diagram gives an 
opposite UV divergence that cancels the IR divergence of the pNRQCD potentials, 
yielding a finite spectrum. 

2. Alternatively one can observe that the pole of the divergence is exactly opposite to 
the infrared pole of the pNRQCD potentials, which can be read from |2^ and the two 
therefore cancel. More precisely, the scaleless, and hence vanishing in dimensional 
regularization, integral of the vacuum part of Eq. (|4.3D , with the octet propagator 
expanded at the cubic order, can be rewritten as the sum of an infrared and an 
ultraviolet divergent integral. The infrared pole cancels with the one in Eq. ( 4.16| ) 



coming from the thermal part, whereas the ultraviolet one cancels the IR divergence 
of the pNRQCD potentials. 

The two interpretation schemes are equivalent and produce at the end a finite spectrum, 
which is the relevant observable. 

The evaluation of r'^{E — ho)^r'^ in ( [4. 161 ), in analogy to what has been performed 
previously in Eqs. ( [4.12| ) and ( [1.13[ ), can be read from [24| 



1 ..„..,.! 1 / N^a^ 



3^i ^ _ ^ I ^''c "s /Ar2 I OAT n \ "s 



r\E - horr'—— = —— -^^ ~ {N^ + 2N,Cf) 



2 



E-hs ^ ' E-hs E-hs\ 8 r ^ " ' mr'^ 

m^ m^ I r \ I h — hs 



+mc - 2C,)'^5\r) + AT - {ylt\] ^^ + 



(4.18) 




Figure 2: Contribution from the NLO chromoelectric correlator: the gluon is longitudinal and the 
dashed blob is the one-loop self-energy. 

where the dots stand for wave function renornializations. Matching to the right-hand side 
of Eq. (|4.9D , we obtain the corresponding contribution to the singlet potential 6Vs of 
pNRQCDhtl: 



Stt V 8 r mr 

m^ m^ I r 



+4(iV, - 2Cp)^6Hr) + iVe^ V^, i [> ) . (4.19) 



Using first-order quantum-mechanical perturbation theory and the value of the Coulomb 
wave function at the origin, IV'n.KO)! = Sio/{t^'>T'^0'o), we obtain the shift of the energy levels 

(cubic) _ Er^Ira^ [ 4g3^ro ^ 8 _ J_ _ 26io\ ^ ^N^^Cp ^ Nl 



''^"'^ " Svr 1 n + ^^^^^ U(2/ + 1) n^ n J ^ n{2l + 1) ^ A 

(4.20) 

Another possible contribution to the potential up to order ma^ is given by radiative 

corrections to the diagram shown in Fig. |l|. At the next order in Os, corresponding to two 

loops, a sizable number of diagrams appears. In [Q, it was shown that in Coulomb gauge 

only one diagram needs to be considered. It consists of a one-loop self-energy insertion 

in the longitudinal part of the chromoelectric correlator and it is shown in Fig. y. It 

contributes at order asTm?j^r'^ , corresponding, in our scale hierarchy, to a magnitude in 

between mg^ and mg^^. Therefore, this term contributes to the spectrum up to order mOg 

only if mg^ ^ T > mg^^'^ . This makes clear that non-static contributions that were not 

considered in the analysis of Ref. Q , such as vertices originating from the spatial center- 

of-mass covariant derivative in the octet sector and higher-order singlet-octet vertices in 

the \/m expansion (see ||3^), contribute to terms smaller than ma^ only. At our accuracy, 

the diagram we are considering can again be written expanding the octet propagator for 

ko r^ T ^ (E — ho) and retaining only the first term, independent of -E — ho- Therefore, 

the result is the same as the one derived in [g]. It reads 

5y^(2 loops) ^ _3^(3) Cp^r'Tml + -C(3) N.Cp al r^ T^ 
2 vr 3 



+i 



^a,r^Tn.^--+,, + lnvr-ln-^ + --41n2-2i^j 



4-7r 



(4.21) 
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where C is the Riemann zeta function (C(2) = '?i"^/6) and the Debye mass rriD is defined as 



mn 



{Nc + Tprif). 



(4.22) 



Equation ( [4.21| ) contains an imaginary part. It comes from the imaginary part of the 
gluon self-energy, which is related to the Landau-damping phenomenon, i.e. the scattering 
of particles carrying momenta of order T in the thermal bath with virtual, space-like 
longitudinal gluons. Furthermore, the imaginary part is infrared divergent. In the EFT 
framework, this divergence has to be cancelled by an opposite ultraviolet divergence coming 
from a lower scale. In the following section, we will indeed show that the same diagram, 
when integrated over momenta of the order of the binding energy, yields the desired UV 



divergence. Finally, we remark that the result in Eq. ( 4.21 ) comes from dimensionally 
regularizing only the integral over k while keeping the thermal part of the gluon self energy, 
which is finite, in exactly four spacetime dimensions. Using the same regularization when 
calculating the contribution coming from the binding-energy scale guarantees that the 
final result for the spectrum is finite and scheme independent. This is not the case for the 
potential, however, whose expression depends on the adopted scheme. 

The contributions to the energy levels and to the thermal width can be obtained 
easily from Eq. (4.21) by using the expectation value for r^ on Coulombic states, i.e. 
(r2)^ ^ = aln^ [Sn^ + 1 - 3/(/ + 1)] /2: 



5E, 

r 



^0' 

(2 loops) 
n,l 

(2 loops) 
n.l 



-^C{3)CF^Tml+^-^N,CFa!T' 
4 vr o 



2 2 

ar,n 



[5n2 + 1 - 3/(/ + 1)] , (4.23) 



Cf 2/^2 

6 ^ 



-+7£; + lnvr-ln^--h--41n2-2 



A'TT 

In 2 iVcCira^T^ 

9 <= -^ >= 



rj.2 2 

/i" 3 

aln^ [5n2 + 1- 3l{l + 1)] 



cm 



(4.24) 



4.1 Summary 

Summing up Eqs. ( 4.14 ), ( 14.19| ) and ( 4.21 ) we obtain the thermal correction to the potential 
in pNRQCDjjTL up to terms whose contribution to the spectrum is smaller than ma^: 



vr 



27r 



Ws = -NcCf aiT^r + ^—Cp Os T^ + 



-.2 , OisCpIx 



3m 



37r 



N^a: 



-^^-{N^ + 2NcCf) 



at 



r 






mr 
1^ ' 



-^Ci3) Cf ^r'T ml + las) N^CFa^r'T' 



+t 



^F 2rr 2 

b 



T2 2 



-- + 7S + Invr - In ^5- -F - - 41n2 - 2 






Att 



In 2 N^Cf at r^ T 



2 „2 rp3 



(4.25) 



where the first two terms come from the linear part of Fig. ||, the terms in square brackets 
come from the cubic term and the last three lines originate from the diagram in Fig. |2|. 
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This correction to the potential can be used for T ^ E, m£) no matter what the relative 
size between E and mo is. 

Analogously, the total contribution to the energy levels coming from the scale T is 



5E. 



(T) 



9 2 6m 



+ 



E^Iral [ 4C3<5,o 



37r 



+ N,C^, 



8 



1 



n^ 



2(5/0 



_| £ — 1 1 L 

n(2Z + l) 4 



+ {-la^)Cp^Trnl + \a2>)N,CFalT^'^[f>n^ + l-2,l{l + l)\. 

(4.26) 

The first and the second lines originate from the diagram in Fig. ||, and correspond to 
the linear and cubic terms in the expansion (4.4). The last line originates from the gluon 
self-energy diagram in Fig. ^, which also gives the full contribution of the scale T to the 
width: 






-,(2 loops) 



(4.27) 



5. Contribution to the spectrum from the scale E 

After having integrated out the temperature in the previous section, many different scales 
{E, rriD, Aqcd, . . .) still remain dynamical in pNRQCDjjxL- ^^ o^^ hierarchy, the binding 
energy is much larger than the Debye mass and Aqcd is smaller than all other scales. Our 
purpose is to compute the correction to the spectrum and the width coming from the scales 
E and rriD- This is achieved by computing loop corrections to the singlet propagator in 
pNRQCDjjrpL. We recall that the gauge sector of pNEQCDj^rp^ coincides with the Hard 
Thermal Loop effective Lagrangian. The longitudinal and transverse gluon propagators in 
Coulomb gauge are given in the Hard Thermal Loop effective theory by |31|^ 



Doo^iko,k) 



, 2 , 2/1 ko ko + k±iri 

k^ + mi, 1 m ; 

^ \ 2k ko-k±i'n 



and 

respectively, where 

Ar,a(^o,A;) = - 



R,A^,_ ,,x I r kikj 



Dii' {ko,k) = {^ij--j^ 



Ar,a(A;o,^) 



(5.1) 



(5.2) 



p 1.2 "^1(^0 /,2 .2,ko,_Jko + k±ir, 



2 VA;2 



■In 



2k^ \kQ — k ziz ir] 



± isgn{kQ) rj 



(5.3) 



^The transverse propagator given there contains a misprint: a factor of po/{2p) should be multiplying 
the logarithm in Eq. (27), as follows from the transverse HTL self-energy given in Eq. (17) of the same 
paper. 



12 



and the upper sign refers to the retarded propagator and the lower sign to the advanced 
one. The "11" component can be obtained from the relation (|2.7]). 

We start by evaluating the diagram shown in Fig. |^, whose general expression is given 
in Eq. ( |4.2| ), but now the longitudinal and transverse gluon propagators are given by Eqs. 
(5.1) and ( [5. 21 ). As we shall see, this is the only diagram we need to consider to get the 
spectrum at order ina^. 

At the energy scale, we have Uq ^ {E — ho) and therefore we have to keep the octet 
propagator unexpanded. However two expansions are still possible. 

1. Since k ^ E -^T, the Bose-Einstein distribution can be expanded in 

I T I k , , 

+ TTT7F + ■ ■ ■ ■ (5.4) 



ek/T _i k 2 12T 

2. Moreover, since k ^ E ^ "mo^ the Hard Thermal Loop propagators can be expanded 
in rin?jj/E'^ <^ 1. At the zeroth order, this corresponds to using the propagators given 
in Eqs. ( |2.3[ ) and ( |2.4| ). Some care is required in the expansion of the transverse 
gluons due to a collinear region, as we shall see later on. 

In the following, we will call 6T,s{E) the contribution of the diagram in Fig. |l] to the 
singlet self energy; the corresponding energy shift and width for the state \n,l) are given 
by 5En,i = {nJ\Re6'Es{En,i)\n,l) and r„,/ = -2{n,l\lm5'Ss{En,i)\n,l) . 

We now proceed to the evaluation of Eq. (^]^) for loop momenta of the order of the 
binding energy, with the HTL propagators defined in Eqs. (^]^) and (^^). We find conve- 
nient to compute separately the contributions coming from the transverse and longitudinal 
gluons. 

5.1 Transverse gluon contribution 

The contribution of transverse gluons to Eq. ([4.2|) is in pNRQCDjj-pj^ 



rv^(trans)('p\ _ ,-„2 ^ ^ ^^»,.4-D f d k l 2 

6^ >{E) - -^g Cf ^-jv t, / (^^^^^T^^^fe^^^o 



+ {l + Mko) ] (Ar(A;o, k) - AaI/cq, k)) 



An{ko,k) + AA{ko,k) 
2 

r' . (5.5) 



We start by evaluating the contribution of the symmetric part, which turns out to be the 
leading one: 

a'r ^-^V4-D f d^k kl (T (E 

x[AR(A:o,A:)-AA(fco,A:)]r% (5.6) 

where we have expanded the Bose-Einstein distribution. The expansion of the HTL prop- 
agators for mo <^ kQ,k needs to be performed with care in the region around the light 
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cone, where the gluon propagator becomes singular. We refer to Appendix ^ for details 
on the expansion and the evaluation of the integral, whose final result reads 



(5.7) 
The suppressed term of order a^r^E^ jT comes from the k/{12T) term in the expansion of 
the thermal distribution, whereas the term of order agTrnj^r'^ / E"^ comes from subleading 
terms in the expansion of the propagator.^ 

We now consider the first term in the square brackets in Eq. ( p.5[ ); it does not depend 
on the Bose-Einstein distribution and, when expanded for kQ,k ^ E ^ rno, gives 

^ rvr\ rC 



where P stands for the principal value prescription. Plugging Eq. (^) back into Eq. ( p^ ) 
yields 

g2 Cp r^H^rV^-^ I -^ — 1- kl \iV-r^ + O iml/E^) 

^ D-1 ^ J (27r)^ E-ho-ko + ir] ^ [ k^ - k^ ^ ^' ' 




r\E-hofr' + 0(a^Emjjr') . (5.9) 



."sC'i? j^^ , 3 j , 2 2 



3 

Summing up Eqs. (5.7) and ( [5.9D we obtain the complete contribution of the transverse 
modes 

^5.(trans)(^) ^ _il^^CFTr\E - Kfr' - i^^^T^E - Kfr' + i'^^^^^^^^^^^ 

X Mn2-ij +C'(as^m|)rV-E^asr2^^/^,Qs^m^r2) . (5.10) 

We remark that the contribution of the transverse modes at the energy scale is imag- 
inary and finite, in contrast with what happens at zero temperature, where it is real and 
UV divergent, the divergence cancelling the infrared divergences appearing in the static, 
1/m and l/m? potentials at the scale 1/r. This is related to the discussion made in 
the previous section regarding the cancellation of the IR divergence in Eqs. ( [4.16| ) and 



(4.20) and can be understood in the following way. For E ^ rriD, the Hard Thermal 
Loop transverse propagator can be expanded for small rriD, giving, at the zeroth order, 
(Ar + Aa)/2 = iP[l/(A:g - k"^)] and (Ar - Aa) = 27rsgn(A;o)5(A:g - k"^). When plugged in 
Eq. ( |5.5D we obtain Eq. ( [4. 3D . Evaluated at the binding energy scale, the vacuum part is 



UV divergent and can be read from |28 



n — 9 r d^h ^2 



D-l ^ J (27r)^ E-ho-kQ + irikl-k'^ + ir] 
-r{E - ho) - + 2 In h7£ - 77 -mvr r', (5.11) 



37r \e II 3 



^This term is of order raoti or bigger only in the very tiny window mg'^ ^ T > mg^'^^'^. For this reason, 



we 



will not include terms of order asTrrij^r jE or smaller obtained from the expansion in rrij^/E' 
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where the logarithm of the energy gives rise to the so-cahed QCD Bethe logarithm in the 
spectrum |28, ^]. On the other hand, the temperature-dependent part gives 



-i-as CpTr'iE - hoYr' - -^-^r\E - hof - + 2 In 



-i7rsgn{E -ho)+jE-^- Invr j r' , (5.12) 

where the term proportional to r*(£' — ho)^r'^ comes from the first term in the expansion 
of the Bose-Einstein distribution and the one proportional to r*(£' — /io)^r* comes instead 
from the second term in that expansion, see (p.4|). In the sum of Eqs. (p. 11) and (5.12|) 



the real parts, divergences included, cancel out and the imaginary parts combine to give 
the two 77i£)-iiidependent terms of Eq. ( p.lCl| ). This shows that the binding energy scale 
contribution produces two opposite UV divergences. In terms of the two interpretation 
schemes discussed in the previous section, we may understand the cancellation of diver- 
gences in two possible ways. In the first way, the vacuum divergence in Eq. (|5.11 ) cancels 



the IR divergences of the potentials, whereas the UV matter divergence in ( |5.12 ) cancels 



the IR matter divergence from the scale T in (4.16). In the second way, we consider the 
real part of the potential in pNRQCDjj-pL as finite, the IR divergences from the scales 
1/r and T cancelling each other, and no UV divergences coming from the energy scale. 



which, as shown by Eq. ( ^.10 ), is indeed the case. We stress that the cancellation of the 



divergences between the vacuum and thermal parts in Eqs. ( [5.11| ) and ( 5.12| ) is due to the 



second term in the low- momentum expansion of the Bose-Einstein distribution, i.e. —1/2, 
which is known in thermal field theory to cause cancellations with the vacuum contribution. 
Finally, we observe that an analogous cancellation is also obtained in the Abelian case Q. 
In order to obtain the contribution to the width from Eq. ( ^.lOf ), we need to evaluate 
r*(£' — ho)'^r^- We proceed as in the previous section and rewrite {E — ho)^ as {E — h^)^ — 
{{E- hs),AV} + Ay2. One then has 



Nj 2 , 2iVeas V 

A S 

4 mr m^ 



r\E - KYr^ = -^at + -^^ + -^ + ... , (5.13) 



where the dots stand for contributions that vanish on the physical state. The width thus 
reads 

' 3 3m 3 mn'^ao 

, 2Enal (4C^ 2 f 8 _ J_ _ 2^\ , 'iN^Cr i^ 

3 1 n ^^\n{2l + l) v? n ) n{2l + l) 4 

- ^ 3 ' ^^[5n2 + l-3/(/ + l)], (5.14) 

where the first line is the contribution from the term proportional to r*(£' — ho^r'^, the 
second line comes from the cubic term and has been obtained using Eqs. ( [4.18D and (|4.2C| ), 
and the third line is the contribution from the last term in the first line of Eq. ( ^.lOj ) . 
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The leading contribution to Eq. ( |5.14 ) is given by the first three terms, which are 



of the same size. The first term comes from the static potential and agrees with the one 
calculated in |^]. The second and third terms come from the kinetic energy; the second 
one agrees with the one calculated in ||^]. This contribution to the thermal decay width 
originates from the possible break up of a quark-antiquark colour-singlet state into an 
unbound quark-antiquark colour-octet state: a process that is kinematically allowed only 
in a medium |^]. Clearly, the singlet to octet break up is a different phenomenon with 
respect to the Landau damping, which, in the previous section, provided another source 
for the in medium thermal width. In the situation E ^ rriD, which is the situation of 
interest for this work, the singlet to octet break up provides the dominant contribution to 
the thermal width. Indeed, comparing the Landau-damping width ( [4.24 ) with the singlet 



to octet break-up width ( ^.14 ) , we see that the latter is larger than the former by a factor 
(mag/m^))^. 



5.2 Longitudinal gluon contribution 

The contribution of the longitudinal gluons to Eq. (|4.2|) is 



5E('-)(i^) = -.,^C.-^/-- / fl ^ ' ,^ , . fc^ 



D-r J {2tt)^ E - ho - ko + i7] 



D^^{ko,k)+D^Q{ko,k) 



+ (\ + nBiko) ) {D^ko, k) - D^oiko, k)) 



2 

rS(5.15) 



where Dqq (k) is the HTL propagator in ( |5.1[ ). The first term in square brackets, i.e. 
(Dqq + Dqq)/2, does not depend on the Bose-Einstein distribution; therefore only the 
expansion in m/j <^ E, corresponding to m/j <C ko,k, is possible. We then have {Dqq + 
Dqq)/2 = i/k"^ + O (rnj^/k^). The first term is the free propagator, which gives a scaleless 
integration, whereas the second one can be shown to contribute at order asEmj^r'^ , which 
is smaller than mag. 

For what concerns the symmetric part of the propagator, i.e. (1/2 + nB{k(j)){D^Q — 
Dqq), it should be noted that the retarded and advanced propagators depend on ko only 
through the HTL self-energy; therefore, imaginary parts in their denominators can enter 



only through the logarithm appearing in Eq. (5.1). Hence, the symmetric propagator is 



non-zero solely in the spacelike /c^ > /cq region, which is related to the Landau-damping 
phenomenon. At leading order in the expansions of the Bose-Einstein distribution and of 
the propagator for mj^/k'^ <C 1, we thus have 

Q + nB(fco)) {Dfoiko, k) - D^h. k)) = ^^^^ (A;^ - kl) + O {ml/k\ Tm%/k') . 

(5.16) 
The first term contributes to the spectrum at order agTrnj^r'^ , so further terms in Eq. 
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( 5.16 ) are not needed (see footnote §). We then have 



S^i^o^s)(E) = g'^Cp 



D-l 
6 



rV'-"" 



d^k 



e 



(27r)^ E-ho-ko + ir] k^ 



^"^^9 {k^ - kl) r^ 



2t: sgn{E - ho) + i [ -+ln^ ^^ + ^iE - --\m: 



r . 



(5.17) 



Equation ( 5.17| ) translates into the following shift of the energy levels 



(5.18) 



For what concerns the width, we observe that the divergence is of ultraviolet origin and 
cancels the one in Eq. ( [4.24 ), yielding a finite width; some care is, however, required in 
the handling of the logarithms of the energy, which give rise to an analogue of the Bethe 
logarithm. We have 



r, 



(long) 
n,l 



^ ' - + In 4- + 7s 



^ ln7rj^[5n2 + l-3/(/ + l)] 



,2 /-»2 „,2 



^ 3 El "•' 



(5.19) 



where Ei = —mCpa^/4 is the energy of the ground state and 



'ni 



K 



Clal J (27r) 



^'' ■(n,/|r|k)|^ln ^^ 



En — k"^ /m 



(5.20) 



(n,Z|r|k) is the matrix element between a (bound) eigenstate \n,l) of hg and a continuum 
eigenstate |k) of ho- This expression can be reduced to a single integral using the techniques 
of [28, 24]. We obtain for a singlet nS state and an octet P wave (the matrix element 
introduces a A/ = 1 selection rule) 



'0 



where 



>r"M 



yi-' 



(z.2 + p2,3 n 



(5.21) 



(5.22) 



The definitions of Y^ , X^ for n = 1,2,3 and pn can be found in |28| and [^], the latter 
reference correcting some misprints in the former. A numerical evaluation of these integrals 
for the three most tightly bound I = states yields: 



I^Q = -0.49673, hfi = 0.64070, /3,o = 1.18970. 



(5.23) 
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5.3 Summary 

In summary, the contribution to the energy levels coming from the binding energy scale is 
entirely due to the longitudinal part of the chromoelectric correlator, 

^^i?=<r^ (5.24) 

which may be read from Eq. ( ^.18| ). The contribution to the decay width coming from the 
binding energy scale is the sum of F^ '^^^'^^ and T^"'^^ : 



NfCpaiT CpUsTEn + -NcCpatT 



2E^ r 4C7|5io 2 ( 8 _ J_ _ ^V 2JVfC^ iV| 

3 \ n '^V"(2Z + 1) n2 n y n(2/ + 1) 4 

\e ;U^ 3 / 

2asCFTm2, C^a^ 



where the first two lines come from the first two in Eq. ( ^.14 ) and the last two from Eq 



( 5.19 ) and from the last term in ( 5.14 ). I^^i is defined in Eq. ( 5.20 ). 



6. Contributions to the spectrum from the scale mj^ 

In our hierarchy of energy scales, the next scale after the binding energy is the Debye mass. 
We thus have to evaluate Eqs. (|5.5| ) and ( 5.15| ) for momenta of the order of mj;). In detail, 



we have two regions to analyze: the first one is /cq ~ .E — ho^ k ~ rriD, corresponding to 
having the octet propagator unexpanded and conversely expanding the HTL propagators 
for kQ ^ k. It can be easily shown that both the transverse and the longitudinal parts 
result in a series of scaleless integrations over k, which vanish in dimensional regularization. 
The second region corresponds to having /cq ~ mo and k ~ m£) : the octet propagator 
then needs to be expanded, whereas the HTL propagators are to be kept in their resummed 
form. The resulting integrals are quite involved, however, by power counting arguments, it 
can be easily seen from Eqs. ( |5.5| ) and ( 5.15| ) that, once the octet propagator is expanded. 



the largest term comes again from the symmetric part of the gluon propagator, due to the 
T/Zcq enhancement factor. The size of this term turns out to be of order asTm^r'^ / E and, 
since we have assumed {ttid/E) <C g, it is beyond ma^ . 

7. Conclusions 

We have computed the heavy quarkonium energy levels and widths in a quark-gluon plasma 
of temperature T such that mas ^ T ^ ma^ ^ tud- Assuming {rriD/E)^ <C g, the 
spectrum is accurate up to order ma^. 

The thermal shift of the energy levels induced by the medium is obtained by summing 
the contribution from the scale T, given in Eq. ( 4.26| ), with the thermal part of the 
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contribution from the energy scale. We remark that the contribution from the energy 



scale, given in Eq. (5.24), is the sum of both vacuum and thermal contributions, which, 
in the transverse sector, cancel. The thermal contribution of the transverse modes can be 
derived from Eq. (5.12). The complete thermal contribution to the spectrum up to order 
mUg reads (we recall that £"„ = —mCpag/{An?) and gq = 2/{mCFas)) 



6e: 



(thermal) 
n,l 



In^Cf al T^^ [3n2 - l{l + 1)] + ^C^ af t^ «^ 



+ 



Sir 



log 



2ttT 



27E 



[ 4C3^ro 



+ N,CJ 



n 



^c^F 



+ 



2N^Cf 
n(2/ + i; 



8 l_ 

n{2l + 1) ~ ^ 



2(5/0 



n 



+ 



+ 



2^„C3a3 

n ^n;t 

Sir 



+ ^[5n2 + l-3/(/ + l) 



27r 



C(3) 



Cf O-s T UIf) 



+ -C(3)7VeCi.a2r3 



(7.1) 



where L^^i is the QCD Bethe logarithm, defined as p8| , 24] 

1 f d^k 



L 



nl 



ClalEn 



(27r)3 



(n, / r k) 



En 



In 



^1 



En — k"^ /m 



(7.2) 



We refer to pq , ^J] for details on the numerical evaluation of this integral. We furthermore 
remark that the thermal contribution to the spectrum is finite, the IR divergence in Eq. 
(4.26) having cancelled against the UV divergence coming from Eq. (|5.12| ). 

The thermal width is obtained by summing the contribution from the scale T, given in 
Eq. (4.27), with the one coming from the energy scale as given in ( ^.25] ), the IR divergence 
in the former cancelling against the UV divergence in the latter. We then have 



-^(thermal) 
■ nJ 



\nICf4t+^^^^{Cf + n^) 



n^ 



+ 



2Enal r4C|Ao 



Cf 
6 



n 



+ N,Cl 



E: 



n(2/ + l) 



n^ 



2(5/0 



n 



_| £ — £_ _| £. 

n(2/ + l) 4 



C'(2)\ , Att 



a,Tml In -^ + 27i5 - 3 - log4 - 2^ + — In 2 N^Cf al T 



J^2 



^2);* 



X alin? 



[hv? + 1 - 3/(/ + 1) 



+ o CfQ^s Tm^F) Oq'^^ In,i 



(7.3) 



where In,i is defined in Eq. ( p. 20 ). We remark that, up to the order considered here, the 
thermal contribution to the spectrum and to the width is independent of the spin. 

Our results are expected to be relevant for the ground states of bottomonium (T(IS') 
and r/b), and to a lesser extent to those of charmonium [J /"ip and rjc), for a certain range 
of temperatures in the quark-gluon plasma for which (^]^) is fulfilled. Let us now try to 
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figure out what our results imply for the electromagnetic decays to lepton pairs or to two 
photons. First of all, the masses of the heavy quarkonium states increase quadratically with 
the temperature at leading order (first line of ( [7.1D ), which would translate into the same 
functional increase in the energy of the outgoing leptons and photons if produced by the 
quarkonium in the plasma. Second, since electromagnetic decays occur at short distances 
(~ 1/m <^ 1/r), the standard NRQCD factorization formulas hold, and, at leading order, 
all the temperature dependence is encoded in the wave function at the origin. The lead- 
ing temperature correction to it comes from first-order quantum-mechanical perturbation 
theory of the first term of ( 4.25| ). The size of this correction is ~ n^T'^ / {rn? a^) ■ Hence, a 



quadratic dependence on the temperature should also be observed in the frequency in which 
leptons or photons are produced by the quarkonium in the plasma. Finally, at leading or- 
der, a decay width linear with temperature is developed (first line of ( |7.3|) ), which implies 
a tendency to decay to the continuum of colour-octet states. Hence, a smaller number of 
vector and pseudoscalar ground states is expected to be in the sample with respect to the 
zero temperature case. 
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A. Details on the evaluation of the transverse HTL contribution 

Our aim is the evaluation of Eq. (^.6|). Owing to the symmetries of the retarded and 
advanced propagators and of the Bose-Einstein distribution we can restrict the integration 



in (5^) to positive values of fco. We then have 

Let us define the quantity A = ko — k. There exist two momentum regions that 



contribute to the integral ( A.l ) for ko r^ k r^ E — hg. We call the first region the off-shell 
region. It is defined by 

X^{E-ho), k^{E-ho), (A.2) 
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i.e. the region where the gluon is far from being on sheh. The second region is cahed the 
collinear region. In this region, we have 

\r^-^^, k^{E-ho). (A.3) 

E -ho 

We observe that the cohinear scale m'j^/{E — ho) has, in our energy scale hierarchy, a 
magnitude in between mg'^ and mg^. It is, therefore, smaller that the Debye mass by a 
factor of itid/E <C 1 and still larger than the non-perturbative magnetic mass, which is of 
order g^T, by a factor T/E » 1. For simplicity, we separate the two regions by a cut-off 
A, such that 

2 

<^-""'»'^»(£^' <*-^) 

We start by analyzing the off-shell region. Here k^ — k"^ = X{2k + A) ^ mj^ and we 
can thus expand the retarded propagator propagator in Eq. (|5.2|) as 



m^ 



* 2 \W 


- ik'o - 


-^^)ivin(|±|g;)) 




{kl- 


-A;2 + i7?)2 



Ar(A:o > 0) - ^2 _ ^2 + i^ + (^2 _ A;2 + i^)2 ^ ^ V (^ " Kf 

(A.5) 

Terms contributing to the real part of this propagator and hence to Ar, — Aa can come 
either from the poles of the denominators, yielding a 8{kQ — k"^), or from the imaginary 
part of the logarithm. However, ^(/cq — /c^) = over the whole off-shell region. We can 
safely discard these terms and obtain 

(AK-A,)(.o>0) = - ">°';;^^°) p^. (A.6) 

Note that the principal value prescription is irrelevant since our integration region excludes 



the poles. From Eq. (A.l), we get 

.^(trans,symm).p. _ g^CpmlT^r^D - 2) [ d^-^k 1 

'^^s.off shell l^J- 2{D-l) J {27r)D-ik^ 

f^^^dko I^^iE-ho + uj^, 
Jo 2-K kl-k-^iE-K^irif-kl ■ ^ -'^ 

This integral does not need to be dimensionally regularized, so we can set D = 4 at this 
point and obtain 

. (trans, symm).. _ g^Cpm^Tr^ f^ dkp 2{E - K + IT]) \ ^ , .^ ( — \ 

'^^.,ofi shell (^)- i27r io 2^ {E-K + tr^Y-kir h^ \kj 

(A.8) 
We consider, now, the collinear region. We start again from the retarded propagator 
introduced in Eq. (|5.2D . We perform the change of variables ko — k = X and we expand for 
A ~ mjj/k <^ k, thereby implementing the collinear hierarchy. We then have 

(Ar - Aa)(A:o > 0) = Ai + A2 + A3 + A4 + A5 + O (^j , (A.9) 
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where the A,- are defined as 



^^ = 2^ 






zrjy 



(A.IO) 



Ao 



64/c4 



,(A 



4fc 



+ ir?)2 (A 



4fc 



i?])^ 



(A.ll) 



A, 



Z771 



D 



In 



2fc 



In 



2fc 
X—irj 



8A;3 



4fc 



+ iry A 






(A.12) 



A4 



im 



In 



D 



2fc 
A+ir? 



In 



2fc 
X—irj 



32k^ 



(A 



^+^^)^ (A 



4fc 



iry)^ 



(A.13) 



A, 



zm 



D 



8A;3 



-W + ^'l A-^-Z7?, 



We start by plugging Ai in Eq. (A.l). We then have 



(A.14) 



^^(trans,symm).^N _ 9 ^F ^j 



d^k 



(2vr) 



A 



w/^c^+^K^-^ 



2{E -ho + ir}) 



{E-ho + iTiY -{k + Xf 



-i-a^CpTr'iE - hofr' + 



a 



.TnA 



D' 



E^ 



(A.15) 



The contribution of A2 is 



<5S 



(trans, symin) 
s,2 



(E) 



g'^CpTr' f d?k f^ dX3im\ 



(27r)3 



'-D 



A 27r 32/c3 
1 



2(^ - /lo + ir/) 



V(A-!j + ir?)2 (A-!J-^r/)2y (ii^ - /^o + ir?)2 - (/c + A)^ 
O {asTm%r^/E^) ; 



the leading order term in the expansion of ((.E — ho + irj)'^ — {k + A)'^)~^, which would 
contribute at order agTrnj-ir^, vanishes because the integral over A is zero. The contribution 
of A3 is 



<5S 



(trans, symm) 
s,3 



(E) 



ig^CpTrnj^r' f d^k 1 



In 



12 

2k 
T 



dX 



2{E ^ho + iri) 



(27r)3 A;2 J ^^ 27r {E-ho + iv? - (^ + A)^ 



-2iTT6 X 



m 



D 

Ak 



2i7r(9(-A)P- 



A 



"-D 
4fc 



r . 
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We then have 



(5S 



(trans, symm) 
s,3 



(E) 



g'^CpTmlr' f^ dk 2{E-ho + ir]) 
127r 7o 2^{E-ho + irjY - k^ 



.n,|| + 



(A.16) 

where the dots mean terms suppressed by 1/A. We now combine this result with the 
contribution from the off-shell region in Eq. (|A.8D to obtain 



i-r-,(trans,syinm) / j;i\ (-r-,(trans,symm) / j-,\ 
OZj„ „„„!,„„ (^£/J + OZj^ 3 [I^) 



s,off shell 






{E-ho + ivy - k^ 
ln2 + ..., (A.17) 



where the dots stand for higher orders. We remark that the dependence on the cut-off 
scale A has disappeared. The contribution of A4 is 



<5S 



(trans, symm) 
sA 



(E) 



ig^Tm^Cpr'^ f d^k 1 



48 
Inl^l 



dX 



2{E -ho + ill) 



(27r)3 fc3 J^ 27r (E - K + irjV - (k + A)^ 



In If I \ 2iTTe{-X) 



(X-^ + zrjy (A 



4k 



17])"^ 



(A - 1^ - ^V)' 



r . 



The needed A integrals are 



and 



^ dX^ 
— m 

-A 2^ 


2k 
X 


/ 1 




1 


[{X-^ + zvr 


(A- 


*-^^)^ 




—i 


l\x \'^ 


ir/)2 


i4:k 

2- + - 



Ak 



m 



2 ' 
D 



so that ST,]. ^^°'^' ^y™™^ ^^~) ]-^g^g only contributions that are suppressed by powers of 1/A. 
Finally, the contribution of A5 is 



6T 



(trans, symm) 
s,5 



(E) 



g^CFr\D-2) f d'^-'k 



2{D-l) J (27r)^-i 



^ dXTiTml 
_a27 2fc2 

X 



^^(A-^) 



mj 

4k 

2{E -ho + iri) 



6 



+ ... 



{E-ho + ir,Y -{k + A)2 

(A.18) 



where the dots stand for higher orders. The contribution of the symmetric part of the 
transverse propagator is then given by the sum of Eqs. ( A.15|) , ( A. 17] ) and (|A.18 ). 
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